Anomalous scaling of passive scalars in rotating flows 
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We present results of direct numerical simulations of passive scalar advection and diffusion in 
turbulent rotating flows. Scaling laws and the development of anisotropy are studied in spectral 
space, and in real space using an axisymmetric decomposition of velocity and passive scalar structure 
functions. The passive scalar is more anisotropic than the velocity field, and its power spectrum 
follows a spectral law consistent with ~ k^^^^. This scaling is explained with phenomenological 
arguments that consider the effect of rotation. Intermittency is characterized using scaling exponents 
and probability density functions of velocity and passive scalar increments. In the presence of 
rotation, intermittency in the velocity field decreases more noticeably than in the passive scalar. The 
scaling exponents show good agreement with Kraichnan's prediction for passive scalar intermittency 
in two-dimensions, after correcting for the observed scaling of the second order exponent. 

PACS numbers: 47.32.Ef; 47.51. +a; 47.27.T-; 47.27.ok 



I. INTRODUCTION 

A passive scalar is a quantity diluted in a fluid in such 
a low concentration that it does not affect the flow evo- 
lution, but that it is advcctcd and diffused by the flow. 
Examples are given by colorant dye used in experiments, 
and aerosols and pollutants in small concentrations in 
the atmosphere. In recent years, the study of passive 
scalars has been associated with substantial advances in 
our theoretical understanding of turbulent flows [H-Q- 
Passive scalars share similarities with three dimensional 
hydrodynamic turbulence, developing a direct cascade (a 
transfer of variance towards smaller scales with constant 
flux), and intermittency (the spontaneous development 
of strong gradients at small scales). As in the study of 
Navier-Stokes turbulence, some topics of interest include 
the persistence of anisotropy at small scales [3, 01 , uni- 
versality, and deviations from scalc-invariance associated 
with intermittent events. For the passive scalar, signifi- 
cant advances to understand these topics have been made 
in the framework of the Kraichnan model [H, 0| • For a 
random delta-correlated in time velocity field, the scaling 
exponents of the passive scalar were obtained and shown 
to have anomalous (intermittent) behavior but to be uni- 
versal. Numerical simulations and experiments showed 
good agreement with these results The results were 
later also extended to consider the behavior of probabil- 
ity density functions, or differences between passive and 
active scalars (see, e.g., [1]). 

Passive scalars in rotating turbulence have received 
less attention. It is well known that one of the effects 
of rotation is to modify the energy transfer, transfer- 
ring the energy preferentiall y to wards modes perpendic- 
ular to the rotation axis |9l4ll|. This results in a two- 
dimensionalization of the flow, and in the development 
of an inverse energy cascade. High resolution simula- 
tions agree with these results, and also indicate that in- 



termittency in the velocity field is substantially reduced 
by rotation [l^ - [T5| . Passive scalar transport has been 
studied in numerical simulations showing that its trans- 
fer is affected by rotation [l^, and in experiments 
[isj showing that passive scalars in rotating flows still 
develop anomalous scaling in their structure functions. 

In this paper we analyze data from direct numerical 
simulations of the Navier-Stokes equation in a rotating 
frame, together with the advection-diffusion equation for 
the passive scalar. Spatial resolution is 512'^ grid points 
in a regular periodic grid. Passive scalar is injected into 
an initially homogeneous and isotropic turbulent flow, 
sustained by an external (random) mechanical force. Af- 
ter reaching a steady state, rotation is turned on. Two 
different scales are considered for the injection of mechan- 
ical energy and passive scalar. In one case, both are in- 
jected at the largest available scale in the domain. In the 
other, both are injected at an intermediate scale, to al- 
low the mechanical energy develop some inverse transfer 
and to see whether this transfer affects the intermittency 
of the passive scalar. Two different rotation rates are 
considered, in both cases chosen to study the regime of 
moderate Rossby numbers. Inertial range scaling is stud- 
ied considering the energy and passive scalar spectra and 
fluxes. Velocity and passive scalar structure functions are 
computed using an axisymmetric decomposition, and the 
corresponding scaling exponents are considered to char- 
acterize intermittency in each field. Finally, probability 
density functions of increments of the velocity and the 
passive scalar are studied, together with visualizations 
of the fields. We find that the passive scalar is more 
anisotropic than the velocity field at small scales, and 



follows a spectral law consistent with ~ k 



-3/2 



where 



denotes wave vectors perpendicular to the rotation 
axis. While intermittency in the velocity field is strongly 
decreased in the rotating case, the passive scalar is still 
intermittent and its scaling exponents show good agree- 
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ment with Kraichnan's prediction for passive scalar in- 
termittency in two-dimensions. 



II. SETUP AND THEORY 

A. Basic equations, code, and simulations 

For an incompressible fluid with uniform mass density 
in a rotating frame, the Navier-Stokes equation for the 
velocity field u, and the equation for the passive scalar 9 
are, 

i9fU + u- Vu = -2r2 X u- Vp + i^V^u + f, (1) 



Vu = 0, 



(2) 



(3) 



where p is the pressure divided by the mass density, v 
is the kinematic viscosity, and k is the scalar diffusivity. 
Here, f is an external force that drives the turbulence, 
(j) is the source of the scalar field, and il = Viz is the 
rotation angular velocity. 

For the analysis in the following sections we use data 
stemming from direct numerical simulations of the above 
equations. We solve Eqs. ([T|), ^ and ^ using a paral- 
lel pseudospectral code in a three dimensional domain of 
size 2tt with periodic boundary conditions [l^ [l^ . The 
pressure is obtained by taking the divergence of Eq. ([T|) , 
using the incompressibility condition ^ , and solving the 
resulting Poisson equation. The equations are evolved in 
time using a second order Runge-Kutta method. The 
code uses the 2/3-rule for dealiasing, and as a result the 
maximum wave number is kmax = -^/3, where TV is the 
number of grid points in each direction. All simulations 
presented are well resolved, in the sense that the dis- 
sipation wave numbers fc^ and are smaller than the 
maximum wave number k^ax at all times. 

The practice of numerically solving flows in a rotating 
frame in periodic boxes dates back to Ref. [2l| . Basically, 
strict periodicity in all three spatial directions replaces 
the hypothesis of homogeneity. With these boundary 
conditions, the centrifugal force (not written in Eq. [T]) 
can be written as a gradient and absorved into the pres- 
sure, and is thus automatically taken into account when 
the Poisson equation is solved. The rotating flow is thus 
considered infinite, homogeneous, with constant rotation 
rate VI, and locally expanded in the rotating frame us- 
ing a Fourier series (in the same way local instabilities 
are studied in unbounded flows) . There are however two 
differences between an infinite and a periodic flow. A 
periodic flow is bounded in the sense that eddies larger 
than the size of the box (27r in our dimensionless units) 
cannot develop. This will be important when we consider 
separation of scales in the inverse energy cascade. Also, 
there is a discussion concerning whether decoupling be- 
tween two-dimensional modes and fast waves takes place 



in rotating flows. While first-order decoupling was shown 
for periodic flows (which have discrete wave numbers) 
[2^ . for infinite domains with continuous wave numbers 
decoupling does not hold [2^. In practice, both idealiza- 
tions have limitations, as other effects in bounded flows 
(e.g., Ekman layers) are not present with these boundary 
conditions. 

The runs arc characterized by Reynolds, Peclet, and 
Rossby numbers. The Reynolds, Schmidt, and Peclet 
numbers are defined as usual as 



Rp — 



UL 



V 

K 



(4) 



(5) 



(6) 



where U is the r.m.s. velocity, and L is the forcing scale 
of the flow defined as 2n/kp, with kp the forcing wave 
number (when the forcing is applied in a wide band of 
wave numbers, fcp- is taken as the minimum of the wave 
numbers in the band). For the simulations, U ~ I, and 
all runs have v = k (i.e., 5'c = 1 and Pg = Re)- 

To characterize the strength of rotation, we use the 
Rossby number 



Rn — 



U 
2Ln' 



(7) 



It is also useful to introduce a micro-Rossby number de- 
fined as the ratio of the r.m.s. vorticity tJrms (w = V x u) 
to the background vorticity. 



Ru. 



2n 



(8) 



For rotation to be important, the Rossby number Ro 
should be smaller than one, but the micro-Rossby num- 
ber should be close to one or larger. If R^ is much 
smaller than one, non-linear interactions are rapidly 
damped by the scrambling effect of Rossby waves, re- 
sulting in a strong quenching of turbulence by rotation 

Several simulations were done with fixed linear spatial 
resolution {N = 512), and same kinematic viscosity {ly = 
K = 6 X 10~^). Parameters for all runs are given in Table 
m The forcing used for the velocity field as well as for the 
passive scalar is a superposition of Fourier modes with 
random phases, delta-correlated in time, and injected at 
the same wave number kp for both fields. One set of runs 
(set A) has forcing applied at fee [1, 2] (therefore kp = 1, 
and the simulations have the largest possible separation 
of scales in the direct energy cascade), while another set 
has forcing a.t k = kp = 3 (set B). This latter choice for 
the injection wave number leaves some room in spectral 
space for an inverse cascade of energy to develop, but also 
reduces the Reynolds number by a factor of three, as scale 
separation between injection and dissipation is smaller. 
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TABLE I: Parameters used in the simulations, kp is the forc- 
ing wave number, Q. the rotation, Ro the Rossby number R^, 
the Reynolds number, and R^j the micro-Rossby number. 



Run 


kp 


n 


Ro 




R^ 


Al 


1 





oo 


1000 


oo 


A2 


1 


4 


0.04 


1000 


0.9 


Bl 


3 





oo 


240 


oo 


B2 


3 


12 


0.04 


240 


0.4 



This results in narrower inertial ranges for all simulations 
in set B. However, the scale separation between the box 
size and the forcing scale will allow us to compare runs 
with and without inverse energy transfer (respectively, 
sets B and A) , and to see whether this transfer has some 
effect in the scaling of the passive scalar. 

The procedure followed in the simulations in both sets 
was the same. A simulation of the Navier-Stokes equation 
with = was done first, until reaching an isotropic and 
homogeneous turbulent steady state (this takes approxi- 
mately ten turnover times). Then, the passive scalar was 
injected, and the run was continued for other ten turnover 
times until reaching a steady state for the passive scalar 
(these runs correspond to run Al of set A, and run Bl of 
set B). Finally, rotation was turned on. Different values 
of O were considered, to have similar Rossby numbers in 
both sets. These runs were started from the last snap- 
shot of the velocity and the passive scalar of run Al or 
Bl. Each of these runs was continued for over twenty 
turnover times. Rossby numbers for each run are listed 
in Table m 



B. Analysis 

Characterization of the flow and passive scalar 
anisotropy, scaling laws, and intermittency is done con- 
sidering power spectra, fluxes, structure functions, and 
probability density functions of field increments. 

Isotropic energy spectrum and power spectrum of the 
scalar field arc defined as usual (summing the power of all 
modes in Fourier space over spherical shells), and denoted 
by E{k) and V{k) respectively. Since rotation introduces 
a preferred direction, anisotropies will be characterized in 
spectral space using the so-called reduced spectra. The 
reduced perpendicular energy spectrum E(k±) and scalar 
power spectrum V{k±) result from summing the power of 
all modes in cylindrical shells of radius k± , with their axis 
aligned with the rotation axis z. The reduced parallel 
spectra -E(fc||) and V{k\^) result from summing the power 
of all modes in planes with k^ = fc|| . Detailed definitions 
can be found in [Tsj . 

Two-dimensional axisymmetric spectra can be also de- 
fined (see e.g., [nj), and give more detailed information 
of spectral anisotropy. Instead, we will consider here 
an axisymmetric decomposition of structure functions, 
which will also give us information of intermittency in 



the velocity and scalar fields. Longitudinal increments of 
the velocity and passive scalar fields are defined as: 

Mx,l) = [u(x + l)-u(x)]-i, (9) 

<5^(x,l) = [0(x-l-l)-0(x)].l (10) 

where the increment 1 can point in any direction. Struc- 
ture functions of order p are then defined as 

Spil) = {duP{^,l)), (11) 

for the velocity field, and as 

Tp(l) = (50''(x,l)), (12) 

for the passive scalar field. Here, brackets denote spacial 
average over all values of x. 

These structure functions depend on the direction of 
the increment. In simulations without rotation, the field 
is isotropic and the 50(3) decomposition was used to cal- 
culate the isotropic component of these structure func- 
tions [24| - |27j . The decomposition was implemented com- 
puting (and averaging) the structure functions for 146 
different directions that cover almost isotropically the 
sphere, using the method described in (28| . In runs with 
rotation, given the preferred direction and the axisym- 
metry associated with it, we will be interested in in- 
crements parallel and perpendicular to fl. We denote 
increments in these two directions as Z|| and l± respec- 
tively. We then follow the procedure explained in [l^ 
to perform a decomposition of structure functions based 
on the 50(2) x K symmetry group (rotations in the xy 
plane plus translations in the z direction). Velocity and 
passive scalar structure functions were computed using 
26 different directions for the increments 1, generated by 
integer multiples of the vectors (1,0,0), (1,1,0), (2,1,0), 
(3,1,0), (0,1,0), (-1,1,0), (1,2,0), (-2,1,0), (-1,2,0), 
(1, 3, 0), (—3, 1, 0), (—1, 3, 0), and (0, 0, 1) for translations 
in z (all vectors are in units of grid points in the simula- 
tions), plus the 13 vectors obtained by multiplying them 
by —1. Once these structure functions were calculated, 
the perpendicular structure functions Sp(l±) and Tp{l±) 
were obtained by averaging over the 24 directions in the 
xy plane, and the parallel structure functions 5p(Z||) and 
Tp{l\\) were computed directly using the generators in the 
z direction. 

For runs in set A, this procedure was applied to eight 
snapshots of the velocity and passive scalar fields sepa- 
rated by at least one turnover time each, while for runs 
in set B we analyzed six snapshots. For large enough 
Reynolds number, the structure functions are expected 
to show inertial range scaling, i.e., we expect that for 
some range of scales Sp ~ l^" and Tp ~ l^" {I may be 
replaced by l±_ or /|| in the rotating case), where ^p and 
Cp are, respectively, the scaling exponents of order p of 
the velocity and scalar fields. Scaling exponents shown in 
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FIG. 1: Energy (solid) and passive scalar (dash-dotted) spec- 
trum for run Al. Kolmogorov scaling is shown as a reference. 



the following sections are calculated for all the snapshots 
analyzed in each simulation, and averaged over time. Er- 
rors are then defined as the mean square error; e.g., for 
the passive scalar exponents, the error is 



N 



(13) 



i=l 



where N is the number of snapshots of the field, d^p- is the 
slope obtained from a least square fit for the i-th snap- 
shot, and is the mean value averaged over all snap- 
shots. The error in the least square calculation of the 
slope for each snapshot is much smaller than this mean 
square error. Extended self-similarity [2^, [s^l is not used 
to obtain the scaling exponents. 



III. NUMERICAL RESULTS 
A. Spectra and fluxes 

In the absence of rotation, velocity and passive scalar 
spectra follow a similar inertial range scaling, close to 
~ /j-s/s expected from Kolmogorov phenomenology, 
except for bottleneck and possible intcrmittency correc- 
tions (the latter only visible in the spectrum at higher 
resolution [Ullll). As an example. Figure [1] shows the 
energy and passive scalar spectra for run Al in the turbu- 
lent steady state. The slope indicates a ~ scaling 
law. 

Once rotation is turned on, the spectral scaling of the 
passive scalar and velocity changes. Figure ^a.) shows 
the reduced perpendicular energy spectrum E{k±) for 
run A2, together with the reduced perpendicular passive 
scalar power spectrum V{k±). The spectral scaling of 
the velocity field changes and gets closer to the expected 



E{k^) ^ 

1.-3/2 



fcj^ [1^ |33[ , while the passive scalar is close to 
scaling. This can be further confirmed in the 
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FIG. 2: (a) Energy (solid) and passive scalar (dash-dotted) 
reduced perpendicular spectrum for run A2. Slopes are shown 
as a reference, (b) Reduced perpendicular energy spectrum 
compensated by kj^, and passive scalar power spectrum com- 
pensated by k^^'^ for the same run. 



compensated spectra shown in [D^b) . Figure [3] shows the 
reduced perpendicular energy and passive scalar spectra 
for run B2; E(kj_) is also close to ~ kj^, and V{k±^) 

3/2 

again is close to k^ . Note run B2 is forced at /ci? = 3, 
and as already mentioned, in the presence of rotation 
there is some room in spectral space for inverse trans- 
fer of energy towards larger scales (indeed, at much later 
times the energy spectrum in this run peaks at k = 1). 
However, the scaling of the passive scalar seems to be 
insensitive to the development of such an inverse energy 
transfer. Since both sets of runs show similar scaling, in 
the following we will concentrate on runs in set A (which 
have a larger scale separation between injection and dis- 
sipation, and therefore a better resolved direct cascade 
range), and only resort to runs in set B for comparisons. 

Unlike the scaling laws observed in the reduced per- 
pendicular spectra of runs with rotation, isotropic and 
reduced parallel spectra show less clear behavior. Figure 
lll^a) shows the isotropic energy and passive scalar spec- 
tra for run A2. The isotropic energy spectrum shows a 
scaling law similar to k~^. However, the isotropic passive 
scalar power spectrum shows a steeper spectrum with no 
clear inertial range, hinting at a more anisotropic spec- 
tral distribution than the energy at scales smaller than 
the forcing scale. In the case of the reduced parallel spec- 



5 



TABLE II: Anisotropy in runs with rotation. -B(fc|| — 0)/E 
is the ratio of energy in all modes with fc|| = to the total 
energy, ^(fcii = 0)/V is the ratio of scalar variance in all 
modes with fcy = to the total scalar variance, tan^(au) is 
the square tangent of the Shebalin angle for the velocity field, 
and tan^(ae) is the square tangent of the Shebalin angle for 
the passive scalar. 

Run £;(fc||=0)/g V(fc||=0)/V tan^(Q^) tan^(a9) 



A2 
B2 



0.5 
0.2 



0.4 
0.1 



13 
14 



20 
50 



tra, no clear power law can be identified for E{k\\) and 
F(fc||); see Fig.g^b). 

Several quantities can be used to measure the devel- 
opment of anisotropics in a rotating flow [ll|, [sj] . As 
an example, the ratio of energy in all modes with fc|| = 
to the total energy, i.e., £'(fc|| = can be used to 

characterize large scale anisotropy [ij]- For a purely 
two-dimensional flow, this ratio is equal to one. For the 
passive scalar, the equivalent quantity V{k^^ = 0)/V can 
also be used. In the simulations without rotation, the 
mean values of E{k^\ =0)/E and V{k^\ = 0)/V are close 
to 0.1. As can be seen in Table HIl these ratios increase 
when rotation is turned on, with £^(fe|| = 0)/E larger 
than — 0)/V in all cases. This suggests that, at 

large scales, the velocity field is more anisotropic than 
the passive scalar field. 

To quantify small-scale anisotropy, the Shebalin angles 
can be used instead. These were originally introduced 
to study anisotropy in magnetohydrodynamic turbulence 
(35I [36j , and for the velocity are given by 



tan^ (ati) = 2 lim 



82(1 J 



= 2 



(14) 



The angle gives a global measure of small-scale 
anisotropy, with a value of tan^ = 2 corresponding to 
an isotropic flow. The square tangent of the Shebalin an- 
gles for the velocity and passive scalar fields, respectively, 
tan^(a„) and tan^(ae), arc given in Table HIl for runs with 
rotation. The values of tan^(Q;e) are larger than the val- 
ues of tan^(a„), indicating that the passive scalar field 
is more anisotropic than the velocity field in the small 
scales. This slower recovery of small-scale isotropy by 
the passive scalar, when compared with the velocity field, 
has already been reported in the case of non-rotating tur- 
bulence, see e.g., [5|, but not in the rotating case to the 
best of our knowledge. 

The growth of energy at fc = 1, and the flux of 
energy towards large scales, associated with the two- 
dimensionalization of the flow resulting from rotation, 
can only be observed in run B2 (see Fig. [S]). As men- 
tioned before, this run has some room in spectral space 
between the largest wave number in the box and the forc- 
ing wave number. Figure [5] shows the energy and passive 
scalar fluxes in runs Bl and B2 (without and with ro- 
tation, respectively) as a comparison. In the absence of 



rotation, the energy flux shows a direct cascade range (a 
range of positive and approximately constant flux), while 
the energy flux is negligible for wave numbers smaller 
than the forcing wave number (fc < fc^^). The passive 
scalar flux shows a similar behavior. In the presence of 
rotation, the energy flux becomes negative for k < fcj?, in- 
dicating energy is transferred towards scales larger than 
the forcing scale, and the flux towards smaller scales de- 
creases. Scale separation is too small to talk of an inverse 
cascade of energy with constant flux, but computational 
limitations in our attempt to have a well resolved direct 
cascade of passive scalar limits our ability to resolve an 
inverse energy cascade. For the passive scalar, no flux to- 
ward larger scales is observed, and the amplitude of the 
positive (direct) flux decreases by a small fraction. 

Based on these results, we can put forward a phe- 
nomenological argument considering the effect of rota- 
tion, to obtain the inertial range spectrum of passive 
scalar observed in the simulations. From Eq. ([3]), we 
can estimate the flux of passive scalar as 



(15) 



where ui^ is the characteristic velocity at the scale l±, 
and 61^ is the characteristic concentration of passive 
scalar at the same scale. Using that the passive scalar 
flux is constant in the inertial range, then scales as 
9f^ ^ al±/ui^. If the energy spectrum in the direct en- 
ergy cascade range follows a power law E{k±) ^ kj^, 
then in the absence of intermittency corrections uij^ ~ 

l]/^, and 



y(fc_L) ~ crfc^ 



3/2 



(16) 



B. Structure functions and scaling exponents 

We now consider isotropic and axisymmetric velocity 
and passive scalar structure functions for the simulations 
described above, using the method explained in Sec. Ill Bl 

Figure shows the result of computing the passive 
scalar structure functions in run Al (isotropic) and A2 
(axisymmetric, only structure functions with l±_ incre- 
ments are shown) , for one instantaneous snapshot of the 
scalar fleld (at t = 7 in run Al, and at i = 18 in run A2). 
The structure functions show a range of scales with ap- 
proximately power law scaling at intermediate scales, and 
at the smallest scales approach the ~ P scaling expected 
for a smooth field in the dissipative range. 

Figure [7] shows a detail for rotating runs A2 and B2 
of the passive scalar second order axisymmetric struc- 
ture functions in the perpendicular direction T2{1±_) (i.e., 
averaged in all directions in the xy plane), and in the par- 
allel direction, T2{1\\). Stronger anisotropy is observed at 
small scales, manifested as a different (steeper) depen- 
dence of 22(^11) with the spatial increment when com- 
pared with T2{1±_), and as a much smaller amplitude of 
T2{1\\) at small scales. 
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FIG. 3: (a) Reduced perpendicular energy (solid) and passive 
scalar (dash-dotted) spectra for run B2. (b) Reduced perpen- 
dicular energy spectrum compensated by kj^ and reduced 
perpendicular passive scalar power spectrum compensated by 



An inertial range with power law-scaling can be identi- 
fied at intermediate scales in T2{1±), but not so clearly in 
T2{1\\) (this is specially true for run B2, where the paral- 
lel structure function shows only the smooth scaling /^). 
This is consistent with results shown in the previous sec- 
tion for the reduced passive scalar spectra. Inertial range 
with power-law scaling was observed for the reduced per- 
pendicular spectrum V{kj_), while no clear scaling could 
be identified for V{k\\). Note also that at scales close to 
the forcing scale or larger, the passive scalar distribution 
always seems to be more isotropic, in agreement with the 
results shown in Table HIl 

The slopes indicated as a reference in Fig.[7]correspond 
to the time average of the second order scaling exponent, 
obtained from a best fit in the inertial range of all struc- 
ture functions at different times available for each run 
(eight snapshots for run A2, and six snapshots for run 
B2). The second order scaling exponents (in the perpen- 
dicular direction) for the passive scalar thus obtained are 
C2 « 0.49 ± 0.01 in run A2, and C2 « 0.50 ± 0.01 in run 
B2. 

The values obtained for C2 are in good agreement with 
the V{k±) ~ kj_^^'^ power law found for the reduced per- 
pendicular passive scalar spectrum, which leads on di- 
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FIG. 4: (a) Isotropic energy (solid) and passive scalar (dash- 
dotted) spectra for run A2. (b) Reduced parallel energy 
(solid) and passive scalar (dash-dotted) spectra for the same 
run. 



mensional grounds to T2{1±) • 

Figure [8] shows the axisymmetric second order struc- 
ture functions for the velocity field, for runs with rota- 
tion at late times. Small-scale anisotropy develops, and 
is more clearly seen in run B2 with f2 = 12. At large 
scales, stronger correlations are observed in the parallel 
direction in this latter case, which is to be expected as 
this simulation has some separation between the largest 
scale available in the box and the forcing scale, and the 
flow becomes quasi-2D at large scales. Comparing these 
structure functions with the ones for the passive scalar in 
Fig. [71 it becomes apparent that at the Rossby number 
considered, the small-scale anisotropy associated with ro- 
tation is stronger for the passive scalar than for the ve- 
locity field. 

Inertial range with power law scaling is also observed 
for velocity structure functions. For runs A2 and B2, 
time averaged second order scaling exponents, corre- 
sponding to scaling 82(1 ±) ~ l±, are indicated as a refer- 
ence by the slopes in Figs.[8ja) and lU^b). Second order 
scaling exponents are ^2 ~ 0.96 ± 0.01 for run A2, and 
^2 ~ 0.98 ± 0.01 for run B2. These values are consistent 
with the inertial range scaling found for the perpendicu- 
lar energy spectrum in these runs, E{k±) - kj^, which 
leads to 6*2 (^_l) ^ IJ^- 

From the curves in Fig. [51 scaling exponents can be 
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FIG. 5: (a) Energy flux n(fc) (solid) and passive scalar flux 
a"(fc) (dash-dotted) for run Bl. (b) Same for run B2. 



computed for higher orders. Velocity and passive scalar 
exponents in the direct cascade range were computed for 
runs in set A up to the seventh order. For runs in set B, 
scaling exponents were only computed up to the sixth or- 
der, as those runs have smaller scale separation between 
forcing and dissipation. 

Figure |9] shows the velocity and scalar exponents for 
runs Al and Bl (no rotation). For the velocity field, 
exponents in both runs are similar. The third order ve- 
locity field exponent is ^3 = 0.99 ± 0.02 for run Al, and 
^3 = 0.98 ± 0.01 for run Bl, in good agreement with the 
value of 1 expected from the 4/5-law for isotropic and ho- 
mogeneous turbulence. The velocity field exponents dis- 
play the well known deviations from the linear (^p = p/3) 
Kolmogorov scaling associated with intermittency. This 
deviation from strict scale invariance is often quantified 
in terms of the intermittency exponent /i = 2^3 — ^gj 
which for these runs is [l = 0.24 ± 0.12 for run Al, and 
/i = 0.25 ± 0.06 for run Bl. The higher order velocity 
exponents obtained from these runs are also consistent 
with previous results for non-rotating turbulence at large 
Reynolds numbers [H [H [H . 

The passive scalar exponents for these two runs also 
display similar values, which are consistent (within error 
bars) with the values of the velocity field exponents for 
orders p ^ 1 and 2, and show larger deviations from the 
Kolmogorov scale invariant prediction for larger values of 
p. This is a well known result which indicates that the 



FIG. 6: (a) Isotropic structure functions for the passive scalar 
up to 8th order in run Al, at t = 7. (b) Axisymmetric (only 
for l± increments) structure functions for the passive scalar 
up to 8th order in run A2, at t = 18. 



passive scalar in isotropic and homogeneous turbulence 
is more intermittent than the velocity field P, . For 
advection and diffusion of a passive scalar in a random, 
delta-correlated in time velocity field in a space with di- 
mensionality d, Kraichnan obtained an expression for the 
passive scalar scaling exponents of all orders as a function 
of the second order exponent ^2 [H, 



V2dC2P+(d-C2)2 + (d-C2) 



(17) 



This prediction is also shown in Fig. [5J with the expo- 
nents computed using the value of C2 obtained from the 
simulations, and for d ~ 3. Deviations from this model, 
which start at p = 4, arc small when compared with the 
deviations of the data from the linear scaling. Similar 
deviations were also already reported in previous numer- 
ical simulations of passive scalar transport in turbulent 
flows [1]. 

In the presence of rotation and for small Rossby num- 
ber, velocity scaling exponents of a strictly self-similar 
flow are expected to follow a linear law — p/2, if the 
energy spectrum scales as ^ kj^. For runs A2 and B2, 
velocity scaling exponents are closer to the p/2 linear 
scaling, and display significantly smaller deviations from 
the straight line for the higher order exponents than in 
the non-rotating case (sec Fig. [TO]). Velocity scaling ex- 
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FIG. 7: Axisymmetric second order structure functions for 
the passive scalar: (a) run A2 with = 4 at f = 18, and (b) 
run B2 with SI = 12 at t = 14. T2{1±) corresponds to dashed 
lines and T2{1\\) to full lines. Slopes indicated as references 
correspond to the time average of the scaling exponents, ob- 
tained from a best fit in the inertial range of the structure 
functions at different times. 



FIG. 8: Axisymmetric second order structure functions for the 
velocity field: (a) run A2 with = 4 at t = 18, and (b) run B2 
with Q — 12 a.tt = 14. 82(1 ±) corresponds to dashed lines and 
S2{1\\) to full lines. Slopes indicated as references correspond 
to the time average of the scaling exponents, obtained from 
a best fit in the inertial range of all structure functions at 
different times. 



ponents are similar for both runs with rotation. The in- 
termittency coefficient in runs A2 and B2 is /i = O.litO.l. 
As reported before [H [H, [H [H, , for large enough 
rotation intermittency is substantially decreased. 

Figure [TU] also shows the scaling exponents for the pas- 
sive scalar in runs with rotation. In the scale invariant 
(non-intermittent ) 

Cp = p/4 hnear law is expected 
for passive scalar exponents, if the second order expo- 
nent is C2 = 1/2. Deviations from this linear scaling are 
clearly visible in all runs (see, e.g. Fig. [TO)) . For runs 
A2 and B2, Kraichnan's linear ansatz (|17p adjusts very 
accurately the numerical data with (2 ~ 0.5 and d = 2. 
The good agreement with the model with d = 2 for runs 
with small Rossby number again points to a strong bi- 
dimensionalization of the passive scalar distribution in 
the presence of rotation. 

Intermittency of the passive scalar decreases with the 
Rossby number, although less than for velocity. The 
intermittency coefficient for the passive scalar, /is = 
2C3 - Ce, is Ms = 0.4 ± 0.1 for runs Al and Bl {fl = 0), 




FIG. 9: Scaling exponents (with error bars) as a function of 
the order p, for the velocity (stars for run Al and diamonds 
for run Bl), and the passive scalar (squares for run Al and tri- 
angles for run Bl). The solid line corresponds to Kolmogorov 
scaling = p/3, and the dashed line corresponds to Kraich- 
nan's linear ansatz with rf = 3 and ("2 = 0.62. 
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FIG. 10: Scaling exponents (with error bars) as a function of 
the order p, for the velocity (stars for run A2 and diamonds for 
run B2) and passive scalar (squares for run A2 and triangles 
for run B2). The solid line corresponds to linear scaling = 
p/2, the dash-dotted line corresponds to linear scaling C,p — 
p/4:, and the dotted and dashed lines correspond respectively 
to the linear ansatz with d — 2 and d = 3 using (2 = 0.5. 
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and Ids = 0.2 ± 0.1 and = 0.22 ± 0.09 for runs A2 and 
B2 respectively (n^O). 



C. Probability density functions 

In this section we consider probability density func- 
tions (PDFs) for longitudinal increments and derivatives 
of the a;-component of the velocity field, as well as for the 
passive scalar. In all cases, curves shown are normalized 
by their variance, and a Gaussian curve with unit vari- 
ance is shown as a reference. In a scale invariant flow the 
PDFs of the velocity and scalar increments are expected 
to be Gaussian and to collapse to a single curve when 
properly normalized. On the other hand, in an inter- 
mittent flow, probability density functions are expected 
to have strong non-Gaussian tails with larger amplitude 
as smaller spatial increments are considered. These tails 
arc associated with the largcr-than-Gaussian probability 
of extreme events to occur, which is the signature of in- 
termittency. 

Figure [Til shows the PDFs of the velocity and the pas- 
sive scalar increments for four different values of the spa- 
tial increment I = 1.6, 0.8, 0.4, 0.2, and 0.1, for run Al 
(fi = 0). As a reference, the forcing scale in this runs 
is « 2tt and the dissipative scale is « 0.05; I = 0.8 and 
0.4 clearly correspond to scales in the inertial range. The 
PDFs of velocity and passive scalar increments for ^ = 1.6 
are close to Gaussian, while for smaller spatial increments 
non-Gaussian tails develop. 

Figure[T2]shows the PDFs of passive scalar and velocity 
increments for rotating run A2. For this run deviations 
from Gaussianity in passive scalar PDFs start at I = 0.4, 
and strong asymmetry is observed for I = 0.4, 0.2 and 
0.1. Velocity increments also deviate from Gaussian dis- 
tribution at I = 0.4. The asymmetry in the tails of the 



FIG. 11: Probability density functions for passive scalar and 
velocity increments in run Al, for five different spatial in- 
crements I = 1.6 (solid), 0.8 (dashed), 0.4 (dashed-dotted) , 
0.2 (dashed-dotted-dotted), and 0.1 (long dashes), (a) Pas- 
sive scalar, and (b) i-component of the velocity field. In both 
cases the dotted curve represents a Gaussian distribution with 
unit variance. Smaller increments have stronger non-Gaussian 
tails. 



TABLE III: Skewness (S) and kurtosis (K) for runs Al and 
A2, for different quantities (for scalar or field increments, I is 
the spatial increment considered). 



Quantity 


I 


S{kl) 


5'(A2) 


^(Al) 


K{A2) 


d^e 




0.04 


-0.15 


12.11 


9.01 


59 


1.6 


-0.06 


0.01 


3.1 


2.8 


59 


0.8 


0.006 


-0.04 


3.4 


2.9 


59 


0.4 


0.005 


-0.2 


4.0 


3.3 


59 


0.2 


0.02 


-0.2 


4.6 


4.2 


59 


0.1 


0.03 


-0.2 


6.0 


5.4 






-0.65 


-0.44 


8.15 


7.96 


5ux 


1.6 


-0.19 


-0.08 


2.8 


2.6 


5ux 


0.8 


-0.25 


-0.08 


3.2 


2.8 


5ux 


0.4 


-0.31 


-0.02 


3.6 


3.1 


5ux 


0.2 


-0.36 


-0.05 


4.3 


3.6 


5Ux 


0.1 


-0.44 


-0.2 


5.3 


4.6 



PDFs of passive scalar increments becomes more impor- 
tant when rotation is present. This is in good agreement 
with results found in the previous section, which indicate 
an increment in small scale anisotropy with rotation. 

In figure [T3l we show the PDFs of velocity and passive 
scalar derivatives (respectively dxU^ and dxO)^ in run Al 
{p. = 0). Deviations from Gaussian statistics are ob- 
served for both derivatives. The PDFs of the derivatives 
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FIG. 12: Probability density functions for passive scalar and 
velocity increments in run A2, for five different spatial incre- 
ments I = 1.6 (solid), 0.8 (dashed), 0.4 (dashed-dotted), 0.2 
(dashed-dotted-dotted) , and 0.1 (long dashes), (a) Passive 
scalar, and (b) x-component of the velocity field. A Gaussian 
curve with unit variance is indicated by the dotted curve. 
Again, as intervals are decreased, curves depart more and 
more from the Gaussian distribution. 
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FIG. 13: Probability density functions for velocity and passive 
scalar derivatives in the x direction for run Al. (a) Passive 
scalar, and (b) a;-component of the velocity field. Dotted 
curves represent a Gaussian distribution with unit variance. 



FIG. 14: Probability density functions for velocity and passive 
scalar derivatives in the x direction for run A2. (a) Passive 
scalar, and (b) a;-component of the velocity field. Dotted 
curves represent a Gaussian distribution with unit variance. 



8x9 and dxUx for run A2 (il ^ 0) shown in Fig. [Ml also 
present strong deviations from Gaussianity and asymme- 
try for dxO. Calculation of the skewness for the PDFs 
of dxd shows a value of 5*61 = 0.04 for run Al, and 
Se = —0.15 for run A2. For velocity gradients, skew- 
ness is Su ~ —0.65 for nm Al, and Su = —0.44 for run 
A2. 

Since differences in the tails of the PDFs in the rotating 
and non-rotating cases may not be easily appreciated in 
the figures, Table [nT] lists the skewness and kurtosis for 
several quantities computed for runs Al and A2. For the 
scalar field, the skewness is close to zero for all quantities 
when = 0, but it becomes negative and with larger 
absolute values in the presence of rotation. The kurtosis 
of all quantities studied decreases in the rotating case. 
For the velocity field, the skewness is negative in all cases, 
but substantially decreases in the presence of rotation. 
As for the passive scalar, kurtosis of quantities associated 
with the velocity field also decreases with rotation. 

Many of the results shown indicate both the passive 
scalar and the velocity field become anisotropic in the 
runs with rotations. To illustrate this. Fig. [15] shows 
an horizontal average of the passive scalar, of the z- 
component of the velocity field, and of the ^-component 
of the vorticity, at late times in run A2. Some correlations 
between the distribution of the three quantities can be 
identified, as the large scale configuration is very similar 
in the three cases. However, vertical vorticity develops 
structures at much smaller scales, as expected. A dc- 
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FIG. 15: Horizontal average of the passive scalar (top), ver- 
tical velocity (middle) and vertical vorticity (bottom) for one 
snapshot of the fields at late times in run A2. 



tailed analysis of spatial correlations between the fields, 
and of spatial (turbulent) diffusion of the passive scalar, 
is left for future work. 



neous turbulence. To this end, simulations with spatial 
resolution of 512'^ grid points were performed in a regular 
periodic grid using a pseudospectral method. Turbulence 
was sustained by an external random mechanical forcing. 
The forcing scale was changed to consider the case where 
energy was injected at the largest available scale in the 
computational domain, as well as the case where a small 
separation of scales between the forcing and the size of 
the box allowed for some inverse transfer of energy to de- 
velop in the rotating case. The main focus of our study 
was in the resulting scaling laws in spectral space, and in 
intermittency as measured from structure functions and 
scaling exponents. 

The most important result of this work is that, in ro- 
tating flows at moderate Rossby number, all scaling ex- 
ponents of the passive scalar except for one can be deter- 
mined with good accuracy from Kraichnan's model [l|, Q 
for a flow with dimensionality d = 2 (at least up to the 
highest order considered here, and within the error bars 
of our study). The value of d can be understood from 
the effect of rotation, which transfers the energy prefer- 
entially towards modes perpendicular to the rotation axis 
[9l-[TT| and makes the flow quasi-two-dimensional. More- 
over, all exponents in Kraichnan's model depend on the 
second order scaling exponent for the passive scalar C,2- 
The value of this exponent can be determined from sim- 
ple dimensional arguments, and we showed the value of 
C2 = 0.5 to be consistent with phenomenology and with 
the numerical results. 

The other results include an analysis of scaling of the 
power spectrum of the passive scalar, which in the ro- 
tating case scales as k'^^'^ . This result seems to be 
independent of whether there is scale separation between 
the forcing scale and the largest scale in the box, or not. 
In the former case, it was also observed that while the en- 
ergy develops an inverse transfer in the presence of rota- 
tion, the passive scalar does not. Measures of anisotropy 
based on the spectrum and on structure functions also 
indicate that the passive scalar is more anisotropic at 
small scales (scales smaller than the forcing scale) than 
the velocity, although at larger scales the passive scalar 
is more isotropic than the velocity field. 

Finally, probability density functions of derivatives and 
longitudinal increments of the velocity field and of the 
passive scalar show a decrease of intermittency in the 
rotating case (as indicated, e.g., by a decrease in their 
kurtosis), more pronounced for the velocity than for the 
passive scalar. This is also in agreement with results 
obtained from the structure functions of the scalar and 
velocity field, as well as from the intermittency exponents 
for each field. 



IV. CONCLUSIONS 

In this paper we analyzed data from direct numerical 
simulations of advection and diffusion of a passive scalar 
in a rotating turbulent flow. The results were compared 
with passive scalar dynamics in isotropic and homoge- 
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